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Calculation of light nucleus reaction cross sections in Geant4 

V. Uzhinsky 

CERN, Geneva, Switzerland and LIT, JINR, Dubna, Russia 

(on behalf of the Geant4 hadronic group) 

Total reaction cross sections of light projectile nucleus ( 2 H, 3 H, 3 He and 4 He) 
interactions with nuclei are calculated using Geant4 models, and compared 
with experimental data. It is shown that the models give various predictions 
at low energies, in the region of the Coulomb barrier. "Shen model" (W.-Q. 
Shen et al., Nucl. Phys. A491 (1989) 130) is identified as an improvement 
over other models. 

Five methods/models are used in the Geant4 toolkit [jQ for calculations of nucleus-nucleus cross 
£N| 1 sections: LHEP originated from Geant3, Kox [2], Sihver [3], Shen [4] and Tripathi [6]. They are im- 

plemented in the corresponding classes: G4HadronCrossSections, G4IonsKoxCrossSection, GJ^IonsSihver 
CrossSection, G4IonsShenCrossSection, G^TripathiCrossSection, GJ^TripathiLightCrossSection. The mod- 
els were proposed many years ago. Several experimental measurements of nucleus-nucleus cross sections 
have been published since that time. Thus there is a new possibility to validate the models. Below a 
set of collected experimental data will be presented with the model calculations performed in the Geant4 
framework. The main attention is paid to light projectile nucleus interactions. 

There is also a new phenomenological model proposed by V.N. Grichine in a recent paper [7]. The 
corresponding class is G^GGNuclNuclCrossSection. 

Usually, nucleus-nucleus interactions are characterized by total, elastic and inelastic cross sections. A 
total reaction cross section is the difference between the total and elastic cross sections. The inelastic 
cross section can include cross sections of various reactions. For example, quasi-elastic scattering cross 
sections of processes where a projectile or a target nucleus is left in the ground state and its collision 
^ | partner becomes into an excited states, a cross section of double quasi-elastic scattering of a process when 

projectile and target nuclei are in excited states after a reaction. The excited states can be states with 
excitations of discrete nuclear levels or states in the continuum region. There can be also charge-exchange 
reactions, fusion reactions at low energies, multi-particle production processes at sufficiently high energies 
\ and so. 

The identification of the reactions is experimentally complicated, especially when there can be unob- 
served neutrons in the final states. Thus the comparison of various experimental data is not a simple 
task. Quite often total reaction cross sections are estimated in the optical model after an analysis of 
elastic scattering data. These estimations will be used in the comparison with model predictions. 

Theoretical models are usually phenomenological ones. The optical model, which can predict various 
cross sections, is applied at low energies (below 1 GeV/u). At high energies the Glauber model [51 1§1 [TU] 
is used. In the last decades, the Glauber model with corrections has been applied for an analysis of data 
at low energies. The application of the model for exotic nuclei ( 6 He, 8 He, 11 Li, 10 B, 16 C and so on) 
interactions with nuclei is an active field now. 
1 Experimental data and Geant4 model calculations are presented for d+A, t+A, 3 Hc + A and 4 He + A 

reactions in figures 1-4 . To summarize, 

• The LHEP model gives energy independent cross sections with a sharp threshold at low energy end. 
It does not provide cross sections for 3 He + A reactions. 

• The Sihver model gives energy independent cross sections with a sharp threshold at low energy 
domain. It covers all considered light projectile nuclei. 

• The Kox model gives meaningful results at all energies. The Shen model is an improved Kox model. 

• The Tripathi Light Ion model suppresses cross sections at low energies for light target materials. It 
does not give results for 3 H + A reactions. Its predictions start to deviate from the experimental 
data for heavy targets and 4 He projectile. 

The Shen and Tripathi model calculations arc compared with the experimental data and with pre- 
dictions of a model proposed by C.-T. Liang et al. [TT] in figures 5-8 . There is also calculations by the 
Grichine's model. To summarize, 
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• The Grichine's model calculations deviate from the data in an irregular manner. It is not clear how 
to improve the model. 

• The Liang et al. parameterization describes the data rather well, better than other approaches. Its 
predictions deviate strongly from other model calculations at E/A < 10 McV/u. At the same time, 
the validity of the Liang model predictions for heavy targets is questionable. 

• The Tripathi model was tuned using mainly the optical model estimations (see Fig. 8, 4 He + Ca, Ni 
cross sections). In general, the estimations are below experimental data. 

As seen, all models give different predictions at low energies. It is unknown if any of the models are 
reliable. All the models are available in the Geant4 toolkit, and can be used. 

Two models, Shen and Tripathi, should be considered for a future development. Theoretical estima- 
tions of the cross sections at low and high energies would be useful for this development. 

Because the upper limit of the Tripathi model is 1 GeV/u, it should be well to use the Shen model 
in various physics lists of the Geant4 in the future Geant4 releases. The Tripathi model may be used for 
low energy applications. 

The author is thankful to M. Kelsey, J. Apostolakis, V. Ivanchenko, T. Koi and A. Ribon for interest 
in this work and important remarks. He is also thankful to Jose Manuel Quesada Molina for a very 
important consideration of the cross section behaviour in the region of the Coulomb barrier. 
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Figure 1: d + A total reaction cross sections. Points with error bars are experimental data [TH [T31 HI]. 
Points without error bars are optical model estimations [15] . 



4 





100 1000 

E/N (MeV) 



10000 



2000 - 



1000 - 




10 100 

E |ab /N (MeV) 



1000 




100 1000 

E/N (MeV) 



10000 



4000 r 



3000 



u 

0) 



2000 



1000 




10000 



E |ab /N (MeV) 



Figure 2: t + A total reaction cross sections. Points are experimental data and optical model estimations 
for d + A reactions (see Fig. 1). 
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Figure 3: 3 He + A total reaction cross sections. Points are experimental data 
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Figure 4: 4 He + A total reaction cross sections. Points with error bars are experimental data [T7J [TBI 
[T51 [T4l I^Ol . Points without error bars are optical model estimations [TS] . 
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Figure 5: d + A total reaction cross sections. Points with error bars are experimental data [HI IT51 HPi] . 
Points without error bars are optical model estimations [15] . 
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Figure 6: t + A total reaction cross sections. Points are experimental data and optical model estimations 
for d + A reactions (see Fig. 1). 
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Figure 7: 3 He 4- A total reaction cross sections. Points are experimental data |16| . 
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Figure 8: 4 He + A total reaction cross sections. Points with error bars are experimental data [T7J [TBI 
[T51 [TH HZ01 HH] . Points without error bars are optical model estimations [TS] . 
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